Time-delay and Doppler tests of the Lorentz symmetry of gravity 
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Modifications to the classic time-delay effect and Doppler shift in General Relativity (GR) are 
studied in the context of the Lorentz- violating Standard-Model Extension (SME). We derive the 
leading Lorentz-violating corrections to the time-delay and Doppler shift signals, for a light ray 
passing near a massive body. It is demonstrated that anisotropic coefficients for Lorentz violation 
control a time-dependent behavior of these signals that is qualitatively different from the conven- 
tional case in GR. Estimates of sensitivities to gravity-sector coefficients in the SME are given for 
current and future experiments, including the recent Cassini solar conjunction experiment. 

PACS numbers: 11.30.Cp, 04.25.Nx 



I. INTRODUCTION 

At the present time. General Relativity (GR) remains 
the best known fundamental theory of gravity, describ- 
ing all known classical gravitational phenomena. Exper- 
iments testing this theory spanning 90 years have failed 
to detect any convincing deviations. Despite its continu- 
ing success, there remains widespread interest in pushing 
the limits of experimental tests of GR in order to find 
possible deviations. This is primarily motivated by the 
consensus that there exists a unified fundamental theory 
that successfully meshes GR with the Standard Model 
of particle physics. Such a theory may produce small 
deviations from GR that could manifest themselves in 
sensitive experiments. 

One promising avenue of exploration involves search- 
ing for violations of the principle of local Lorentz sym- 
metry P, J a foundation of GR. Candidate theories ex- 
ist in which this symmetry principle may be broken, at 
least at observable energy scales. These scenarios in- 
clude strings [1, noncommutative field theories Q, 
spacetime- varying fields [6^, quantum gravity 0], su- 
persymmetric theories [8j], random-dynamics models @, 
multiverses [l^ , and brane- world scenarios [llj . 

A general theoretical framework for testing Lorentz 
symmetry in both gravitational and nongravitational sce- 
narios has been developed and is called the Standard- 
Model Extension (SME) The SME is an effective 
field theory that incorporates the known physics of the 
Standard Model and GR, while also including all possi- 
ble Lorentz-violating terms The Lorentz-violating 
terms are constructed from Standard Model and gravita- 
tional fields and coefficients for Lorentz violation, which 
control the degree of the symmetry breaking. 

One useful subset of the SME, called the minimal SME, 
contains the Lorentz-violating terms that dominate at 
low energies. The matter sector of the minimal SME 
has been explored in experimental studies involving light 
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[T 6|. [TtI . [l^ [l9|. electrons I20II . protons and neu- 
trons [21 1 ■ mesons[23| . muons [23], neutrinos [24|, and 
the Higgs [2^. Some nonminimal SME terms, includ- 
ing Lorentz-violating operators of higher mass dimension, 
have already been explored in the photon sector in Refs. 
[2^. In addition, because of the similarities of space- 
time torsion to certain types of Lorentz violation, exper- 
imental searches for SME coefficients have been used to 
place new torsion constraints j2^ . A summary of the cur- 
rent experimental constraints on SME coefficients can be 
found in Ref. 

Studies of the curved spacetime generalization of the 
SME have recently begun. Within the setting of a gen- 
eral Riemann-Cartan spacetime, the dominant SME la- 
grangian terms in the matter and gravitational sector 
have been established • The matter sector of the SME 
couples to gravity via the spin connection and vierbein. 
In this scenario, some novel effects can occur that are 
controlled by certain matter sector coefficients which are 
unobservable in the flat spacetime limit ^29.] . In the pure- 
gravity sector, key experimental signals in the Riemann 
spacetime limit have been established [30l| • Experimental 
work constraining SME coefficients in the gravity sector 
has already begun with atom-interferometric gravimeters 
fs^ . lunar laser ranging Gravity Probe B [s^l, and 
short-range gravity tests [34] . 

Of the classic tests of GR, the so-called fourth test, 
involving the measurement of the Shapiro time delay 
of light passing near a massive body 36], has recently 
gained attention. Improvements in two-way radio com- 
munication with deep-space satellites, such as the Cassini 
probe, make possible a reduction in solar corona noise, 
yielding significant improvements in the accuracy of such 
tests (38j |. Further improvement in both time-delay 
and light-bendiiig tests is also expected in the future 
[H, Ea EH, El S H Si- it is therefore relevant to 
analyze in some detail the signals for Lorentz violation 
in such experiments. Some preliminary results describing 
the leading Lorentz-violating corrections to the Shapiro 
time-delay effect were obtained in Ref. (soj and were ap- 
plied to the case of binary-pulsar timing experiments. We 
seek here to elaborate on these results, determine in ad- 
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dition the associated gravitational frequency shift signal, 
and study potential signals in solar-system experiments. 

This paper is organized as follows. In Sec. [TTl we dis- 
cuss the theoretical foundations of this work. Section 
III Al reviews key results in the gravitational sector of the 
SME, including the post-newtonian metric. We discuss 
light propagation in a background spacetime in Sec. Ill Bl 
and apply the results to obtain the time-delay formula 
in Sec. Ill CI and the frequency shift formula in Sec. IIIDI 
In Sec. mil we examine the results in the solar-system 
scenario. Some preliminary discussion of the experimen- 
tal scenario in Sec. IIII Al ls followed in Sec. IIII Bl bv some 
exploration of the features of the Lorentz-violating sig- 
nals in time-delay tests and Doppler tests. We discuss 
how analysis might proceed and estimate sensitivities for 
existing and future experiments in Sec. IIII CI The main 
results of this work are summarized in Sec. lIVI Through- 
out this work we adopt standard notation and conven- 
tions for the SME, as contained in Refs. [H, [13, BiSI . In 
particular, we work in natural units where h = c = 1 and 
with the metric signature — h -f-t-. 



II. THEORY 
A. Basics 

The SME with gravitational and nongravitational cou- 
plings was presented in Ref. [l^. The general scenario 
is a Riemann-Cartan spacetime and includes couplings 
to curvature and torsion degrees of freedom. We focus 
here on the pure-gravity sector in the Riemann spacetime 
limit, within the minimal SME case. The relevant action 
for this sector of the SME is written as 
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u)R + s'"'Rl 



(1) 



In this expression, g is the determinant of the spacetime 
metric g^j^, R is the Ricci scalar, R^^ is the trace-free 
Ricci tensor, C^x^iv is the Weyl conformal tensor, and G 
is Newton's gravitational constant. The 20 coefficients 
for Lorentz violation u, s^^ , and t'^'^'^^control the leading 
Lorentz-violating gravitational couplings. The additional 
piece of the action denoted S' contains the matter sector 
and possible dynamical terms governing the 20 coeffi- 
cients. 

In the SME formalism, the action maintains general 
coordinate invariance, or observer diffeomorphism sym- 
metry, as well as observer local Lorentz symmetry. How- 
ever, because of the transformation properties of the co- 
efhcients for Lorentz violation, the SME action breaks 
both particle local Lorentz symmetry and particle diffeo- 
morphism symmetry fl^, [Sa ] . In the present context of 
the action in Eq. Ilj, the degree to which the particle 
symmetries are broken is controlled by the coefRcients u, 
st"", and t^^t"". 



It has been demonstrated that explicit breaking of lo- 
cal Lorentz and diffeomorphism symmetry generally con- 
flicts with the Bianchi identities of Riemann geometry 
[l3'|. In the action ([1]) above, explicit symmetry break- 
ing would correspond to specifying a priori the functional 
forms of the coefficients u, s^'', and t'^^^" . If the Lorentz- 
symnietry breaking is dynamical, however, the conflict 
with Riemann geometry is avoided p^ . In the latter 
scenario the coefRcients for Lorentz violation are dynam- 
ical fields and satisfy their own equations of motion. This 
ensures that the Bianchi identities hold. 

We consider here the case of spontaneous Lorentz vio- 
lation. The dynamics governing the coefficients appear- 
ing in Eq. ([1]) are contained in the S' term. Through a 
dynamical process, the coefficient fields acquire vacuum 
expectation values that are denoted as u, s'"', and J'^^^" ^ 
For example, this may occur through the introduction of 
potential terms in S' for u.s^'^, and t'^^^'^ , whose min- 



ima are nonzero [1, U, l35|, |48| . This scenario has been 
treated for the action in Eq. ([!]) in the linearized gravity 
limit, along with a broad study of signals for Lorentz vi- 
olation in gravitational experiments, in Ref. [30.] . In par- 
ticular, the post-newtonian metric was obtained, which 
comprises the starting point of this work. Note that mod- 
els of spontaneous Lorentz-symmetry breaking, capable 
of producing the effective coefficients for Lorentz viola- 
tion in ([1]), exist in the literature. These include scalar 
4i, vector 0, [H, [H [S [s^, [Ml , and two-tensor models 

To study the propagation of light signals in a weak-field 
gravitational system, such as the solar system, the dom- 
inant 0(2) contributions to the post-newtonian metric 
are needed (54] . The relevant terms in the metric for the 
pure-gravity sector of the minimal SME are controlled 
by the coefficients s'^". They can be written in compo- 
nent form, in an asyrnptotically inertial post-newtonian 
coordinate system (Sa), as 

500 = -l + (2 + 3s"°)C/ + sJ'*=C/-'''= + 0(3), 
goj = {ai-2)^^U~ais°''W + 0{3), 

= P'' + [2 + (1 - 2a2)s°"]S^''U + 2{a2 - l)s^''U 

(2) 

In the limit of vanishing s^'^ coefficients, the post- 
newtonian metric of GR is recovered. The potentials 
appearing in this metric are given for an arbitrary mass 
density p by 

U = GI ^^d^x', 



U3k ^ Q 



p{x\t){x-xy{x-x')\ ^ , 



In Eqs. 0, some coordinate gauge freedom remains in 
the two quantities Oi and 02- For example, the standard 
harmonic gauge can be obtained by setting ai = a2 = 1. 
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We leave these quantities unspecified to explicitly display 
the gauge-dependent nature of some of the results we 
derive in this work. As discussed in detail elsewhere [13] , 
the relationship between this metric and the standard 
Parametrized Post-Newtonian (PPN) metric ^] is one 
of partial overlap in a special isotropic limit of the SME. 

We will consider in this work the post-newtonian met- 
ric that is produced by a massive body at rest at the 
origin of the chosen coordinate system. The dominant 
contributions to the potentials appearing in ^ are from 
the monopole terms. They depend on the coordinate po- 
sition of the test body relative to the origin . Thus we 
use 



U 



GM 
r 

GMrh^ 



(4) 



where M is the suitably defined mass of the central body. 
In 11]), we have neglected higher multipoles, which can 
play a role in systematics |45l l46j . and would be needed 
for a full treatment of the general relativistic time-delay 
and Doppler shift signals. For the present purposes, how- 
ever, we need only the dominant contributions to these 
signals that are controlled by the s^"^ coefficients. 

If the mass of the central body is distributed signifi- 
cantly outwards from its center, then a substantial spher- 
ical moment of inertia can arise, as happens with the 
Earth. In this case, for a light signal grazing the surface 
of the central body, terms in the metric proportional to 
the moment of inertia / of the body, as well those that 
might be produced from a quadrupole moment, can give 
a significant contribution to resulting signals controlled 
by the coefficients :30]. For simplicity in this work, 
we neglect such cases and discard the metric terms de- 
pendent on /. This is not expected to produce a severe 
problem in typical solar-system experiments since it is 
known, in terms of the Sun's mass M and radius i?Q, 
that Iq « 0.059Mi?2, [si. 
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Note that under these assumptions we are neglecting 
Lorentz violation in the photon sector of the SME, which 
in any case is tightly constrained compared to the gravi- 
tational sector [il. 111, [13, m, [H . 

We first consider a one-way light signal sent from an 
event E to an event P, that passes a central body. 
We will need to find the deviation of the light ray 
path in curved spacetime from the straight line path in 
Minkowski spacetime between the two events. The spa- 
tial endpoints of the path will be fixed at the two events 
E and P, which amounts to solving ([S|) as a boundary- 
value problem rather than an initial-value problem. 

To find the corrections due to curved spacetime we 
adopt a perturbative method using the linearized expan- 
sion for the metric and an expansion for the wave vector 



(7) 



Here /i^^ are the metric fiuctuations, representing the de- 
viation of g^i, from the flat spacetime metric 77^^. The 
first term in the second equation is the zeroth-order 
wave vector that is constant and satisfies the condition 
rjfj.up'^p'^ = 0. The second term Sp^^ is the correction to 
the wave vector due to curved spacetime. Applying the 
null condition for the full wave vector p^ to leading order 
in the metric perturbation /i^^ yields a constraint on Sp^^: 



-hf^i^p'^p" 



(8) 



We shall denote the coordinates of the endpoint events 
E and P as {te, r^) and (fp, rp, respectively. Generally in 
what follows, quantities referred to each of the events are 
denoted with subscripts e and p. The zeroth-order spatial 
trajectory for the light ray will be a straight line in the 
direction R — rp — f^. This implies that the zeroth-order 
wave vector, tangent to this straight line, has components 
p° = 1 and p> = W, where R = R/ R and R ^ \R\. The 
zeroth-order spatial trajectory can be written as 



4(A) = R^\ + V, 



(9) 



B. Light propagation 

To find both the time-delay signal and the Doppler 
shift signal we employ standard methods and adopt 
the geometric optics limit of electrodynamics in curved 
spacetime [13, [1^. We take the wave vector of a light 
ray, tangent to the light path a:^(A), to be 



p- 



(5) 



where A is an affine parameter. Since the light ray is a 
null geodesic, it obeys the geodesic equation and the null 
vector condition given by 



dpi" 



(6) 



where V is the impact parameter vector. It can be writ- 
ten as 



V = 



Wrr, ■ R. 



(10) 



Furthermore, to be consistent with the boundary condi- 
tions, the parameter A is taken to vary from —le = fe'R 
to Ip — fp- R, from which it follows that le + lp = R- The 
various quantities that we use to describe the zeroth- 
order trajectory of a light ray passing a central body are 
depicted in Fig. [TJ 

The parametrization and definitions above have an im- 
mediate consequence on 6p^ . Integration of the spatial 
components of the definition ^ over the light path, fol- 
lowed by use of the second equation in ([7]) yields 



Sp^dX = 0. 



(11) 
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FIG. 1: Diagram illustrating the meaning of the position vec- 
tors fp, fe, R, and the impact parameter b, for the zeroth- 
order light trajectory passing a massive body from events E 
to P. 



This result is just a reflection of the fact that the spatial 
endpoints of the trajectory are fixed. Equation pTjl will 
be useful in deriving the light travel time formula and 
the Doppler shift formula, as we show below. 

The corresponding integral involving the time compo- 
nent Sp^ does not vanish, however, on account of its being 
fixed by Eq. ([8]). In fact, it can be used to derive the light 
travel time. Integrating the time component of the defi- 
nition ([5]) from E to P, and making use of ([5]) and pT|) . 
we obtain 



te^R 



(12) 



This equation forms the starting point for the derivation 
of the time-delay formula in Sec. Ill CI 

We now consider the shift in the frequency of light 
measured by two observers at the two events E and P. 
The ratio of the frequencies v measured at the two events 
can be obtained from the standard formula 



VE {UePpi)E' 



(13) 



where and are the four- velocities of two distinct 
observers present at events P and E, respectively. Note 
that the quantities in the numerator and denominator are 
to be evaluated at the two events P and E, as indicated. 

To obtain an explicit expression for the frequency shift 
for the one-way trip past a massive body, it will be conve- 
nient to work with the covariant components of the wave 
vector p^ — fffj.^p'^. Expanding Eq. into space and 
time components yields 



dt 
dTr, 



dTe 



Po 



Pq + v^P-i 



(14) 



where and are the coordinate velocities of the two 
observers at events E and P, and Te and Tp are their 
proper times, respectively. 

One convenient consequence of using the covariant 
wave vector is that, for a spacetime metric with no 
explicit time dependence, the component po will be 
constant along the path of light x^{X). Furthermore, 
this condition will be approximately true for the post- 
newtonian metric ([2]) from a massive body approximately 
at rest at the origin of the chosen post-newtonian coor- 
dinate system. Therefore we take 



dpo 
d\ 



0. 



(15) 



It will be important, however, to determine what the 
constant po is, in order to obtain the correct frequency 
shift result. 

To determine the covariant components po and pj of 
the wave vector in Eq. (I14|) we first expand in the manner 
of©: 



Pf^=Pt, + SPf, 



(16) 



where p^ = rj^^p'^ . Using the null constraint Eq. ([5]), 
and the properties oip^, we can establish that 



8po = -W5p^ + hot^p^ - -ht.uP^p" , 



Spj = 5p> + hj^pf". 



(17) 



If we integrate the constant pq = —1-1- Spo over the 
light path, use the first equation in (fTT]). and Eq. (fTTj) . 
we can establish that 



^Po 



1 



1 



R 1 1 ^'^o,r ~ ^h,.rr ] dx. (is) 



Furthermore, if we insert the expansion (|16p into the 
geodesic equation and integrate over the path we 
find 

Sp, (P) - Sp, (E) = i l'"^ djh^^rrdX. (19) 

To find the value of Spj at the endpoints, which is 
needed to evaluate the frequency shift , we start with 
the expression (fTTj) . A suitable integration by parts, fol- 
lowed by the use of (fT9|) . yields the values of Spj at the 
two events P and E in terms of integrals of metric com- 
ponents: 



spae) - ^ / [(^ - ip)djh,,rr + 'ih,^r)]dx. 

(20) 

The expressions ([20l) and ([T8|) form the starting point 
of the derivation of the Doppler shift in Sec. Ill Dl Note 
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that, although we will focus in the next sections on the 
metric from the gravity sector of the minimal SME, the 
results ([12]), and JSO]) could be apphed to al- 

ternative theories of gravity in the linearized limit, with 
an approximately time-independent metric. In particu- 
lar, though it lies beyond the scope of the present work, 
it would be of interest to investigate effects outside of 
the gravity sector of the minimal SME, such as matter- 
gravity couplings [29]. Finally, we note in passing that 
our results in Eqs. (PU|) are consistent with Ref. [58l |. 



C. Time-delay formula 

To establish the one-way light travel time, which con- 
tains a time-delay term due to curved spacetime, we must 
evaluate the integral in Eq. p^ . The projection of the 
metric along that appears in the integrand can be 
written as 



The quantities A and A-''^ are given by 
A 



(21) 



4 + s"°(4 -2a2) + (2ai - A)s°' 
+2{a2 - l)P''Rm\ 

~a2-s^^R^W + 2a2R^R''s°°. (22) 



With these definitions the light travel time takes the form 



if} if 



Ia 
2 



UdX 



2 



W'^dX. (23) 



Using the monopole expressions in Eq. and eval- 
uating the potentials with the zeroth-order spatial tra- 
jectory ([9]), these integrals can be evaluated by standard 
methods. The resulting expression for the one-way light 
travel time, to post-newtonian order 0(2), is given by 



s°^Rn\n 



= R + 2GM{1 
+GM 
+ (02 

+GM[ais°^b> + (as - 2W^R^h 
+..., 



R 



R 



Tp) 



(24) 



where the ellipses represent higher order post-newtonian 
corrections. Neglecting these term suffices to establish 
the leading effects from Lorentz violation for experi- 
ments. Note that this one-way result matches that ob- 
tained in Ref. [sot in the appropriate limit. Also, in the 
isotropic limit of the SME, and for the appropriate coor- 
dinate choice, the result ([M)) matches the standard PPN 
result in, [53. 



In many practical cases, the light signal is reflected 
from a planet or spacecraft. Using we can establish 
the round-trip light travel time. This involves adding the 
light travel time for a signal transmitted by an observer 
at event P that travels to the other observer arriving at 
an event denoted E' . The light travel time for the return 
trip can be obtained from ([24]) with the substitutions 



' pi 

re , 



(25) 



where re ' is the position of event E' . Note that the quan- 
tities R and h will change for the return trip accordingly. 
We assume that the observer at event E, later receiving 
the returned signal at event E' , is traveling at small veloc- 
ities compared to 1. Thus it suffices to approximate the 
motion during the light transit as rectilinear. The small 
velocities are in any case implied by the post-newtonian 
expansion adopted here. If we account for this motion 
during the light transit, but we neglect terms of order 
GMv, the order GM portion of the light travel time is 
equal to its value for the outgoing trip, except for sign 
changes in the s^"* terms. Thus we obtain for the round- 
trip light travel time 



At 



2i?(l -R-v) 



+2GM[-a2S° 



+ 4GAf(l + s™)ln 



R 



R 



+ {a2 - l)P'R^R''] (^-^+^Jl) 
\re rpj 

+2GM{a2 ~ 2)-s^''Wh^ ("".l^lA . (26) 



Note that the terms with the coefficients canceled 
when adding the outgoing and return trip contributions. 
This is due to their oddness under parity. 

Neglecting terms of order GMv, the measured elapsed 
proper time Are at the receiver is related to the above 
result by the factor dTe/dt, which is to be evaluated along 
the worldline of the receiver. In principle, this factor con- 
tains contributions from the spacetime metric near the 
observer present at event E and is related to the classic 
gravitational redshift as discussed in the next subsection. 
For our analysis in this work, we focus on effects from a 
single body stemming from the 0{GM) terms in the ex- 
pression above, though the results could be generalized 
to N bodies. 

There are two key time scales which could be used to 
distinguish the large special-relativistic effects contained 
in the first term in (j26p from the smaller terms of order 
GM jsl]. The time scale over which significant changes 
occur with the first term is essentially an orbital time 
scale r/TJ, where r and v are typical orbital distances and 
velocities, respectively, comparable to R and v defined 
in Sec. Ill Bl The conjunction time scale h/v is approx- 
imately the time scale over which the 0{GM) terms in 
(I26p vary significantly. For typical experiments this is on 
the order of days. 
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The dominant contribution from the 0{GM) terms 
comes from the logarithmic term in (|26p . Note that the 
only coefficient for Lorentz violation appearing in front 
of the logarithmic term is the rotational scalar s^" , which 
points to the possibility of its being measured at the same 
level as the PPN parameter 7. Anisotropic coefficients 
control many of the remaining terms. As we show for 
specific experiments in Sec. IIIIl the typical size of the re- 
maining terms are somewhat suppressed relative to the 
logarithmic term. 

It is important to note that, in principle, the special- 
relativistic terms in (pS)) also receive corrections due 
to the gravity-sector coefficients s'^". These corrections 
would arise through modifications to the orbital dynam- 
ics of the transmitting and reflecting bodies (e.g.. Earth 
and spacecraft or planet). For the purposes of detailed 
modeling, these effects could be included, for example, 
by modeling the orbits as oscillating ellipses. Secular 
changes in the orbital elements due to the coefficients 
s^" could be included using the results from Ref. [30| . 
In any case, such orbital corrections are expected to be 
relevant over the orbital time scale r/v. 



D. Frequency shift 



The result is expanded to all orders in velocity for the 
special-relativistic terms, but to 0(3) in the terms de- 
pending on the metric fluctuations h^^. With some ma- 
nipulation we obtain 



VE 



1 - W2 1 - yOR] 




(28) 



where the term arising from the effects of gravity via the 
metric fluctuations is labeled g and is given by 



g y^Ej us 



vp 

VEj u 



(29) 



The term labeled RS on the right-hand side of ([29]) 
is the gravitational redshift. This term arises from the 
spacetime metric being evaluated at the endpoints of the 
light trajectory, namely events E and P. It is given by 



vp 



= 1 



1 



-00 



GM 



re - rp 



' e' p 



e e 



., (30) 



In GR, in addition to the bending of light and the 
time-delay effect, the frequency of light also changes af- 
ter having passed near a massive body [sj. This effect, 
closely related to the time-delay effect, is distinct from 
the classic gravitational redshift and vanishes for station- 
ary observers. In this section, we evaluate the one-way 
frequency shift, using the results of Sec. Ill Bl and also 
determine the fractional frequency shift for a two-way 
reflected signal. 

We begin with Eq. p4|) . expanded to leading order in 
the wave vector shift 



Up dt dTe 11 — Spo — W — 5pj{P) 
VE dTp dt \l~Spa~ v^W - vUpj{E) 



(27) 



where the ellipses represents higher post-newtonian 
corrections. Equation ([50)1 includes leading Lorentz- 
violating corrections to the standard gravitational red- 
shift of GR, which is recovered in the limit s^'^ — 0. The 
result ((30l) would be of interest to investigate for gravita- 
tional redshift experiments, such as those incorporating 
sensitive atomic clocks on Earth or aboard orbiting satel- 
lites [53,, ,59,, 60, , .62.] . Our main focus in this work, how- 
ever, will be on the time-delay effect and its associated 
contribution to the frequency shift derived below. 

The term labeled D in Eq. ([^5]) is the gravitational 
frequency shift of light due to the wave vector corrections 
Sp^, which reads 



— ) = Spo{v~wyR^ ~w^Spj{P)+v^Spj{E) 
veJ D 



2R 



(31) 



This term represents a gravitational correction to the 
usual Doppler shift of Special Relativity. The integrals 
in (|3ip can be evaluated by inserting the post-newtonian 
metric ^ and using the zeroth-order spatial trajectory 
of the light ray, in a manner similar to Sec. Ill CI The 
result is signiflcantly more cumbersome than (|24|) . and 



so we adopt an approximation that is suitable for cap- 
turing the dominant terms that are proportional to the 
coefficients for Lorentz violation s'^". After evaluating 
the integrals in (pi]) . the results can be grouped accord- 
ing to powers of {GMv /b){h/r)^ . We will focus here on 
near-conjunction time scales where the dominant terms 
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in (|3ip are of order GMv/b and higher order terms will 
be suppressed by powers of the small factor b/r. 

Keeping only the order GMv/b terms, we obtain for 
the frequency shift contribution (|3ip . 



Vp_ 

VE 



AGM , 



[{l + s°° - s"^R^ + s^''Vb^)b - s^^Vb^] . 

(32) 

In this expression the dot denotes a time derivative with 
respect to the post-newtonian coordinate time t. Note 
that the arbitrary quantities ai and 02 keeping track of 
the coordinate gauge freedom have vanished in this re- 
sult, indicating the coordinate invariance of (I32p . The 
result (I32p can also be verified by taking the coordinate 
time derivative of (j24p and using a known relationship 
between the frequency shift and light travel time ^37]- 

The result ([32| can be contrasted with the contribu- 
tions to the frequency shift contained in the remaining 
terms in Eq. and also ([50]) . In the same manner as 
the special-relativistic terms in the light travel time ex- 
pression the velocity contributions in l|28p and the 
gravitational terms in (j30p will vary over the orbital time 
scale r/v in typical experiments. In contrast, the signal in 
P2p will vary most significantly when the light ray passes 
near the central body (6 << r), when the observers and 
the central body are in conjunction. 

We now calculate the fractional frequency shift of a 
light signal reflected from the planet or spacecraft. Thus 
we seek 



Sv_ 

V 



(33) 



where v is the transmitted frequency and v' is the re- 
turned frequency. In a manner similar to what was done 
for the round-trip light travel time in Sec. Ill CI we can 
obtain the frequency shift for the return signal with suit- 
able substitutions in the one-way result ([5^ . Adding the 
return signal to the outgoing one, we find that the grav- 
itational portion of the leading fractional frequency shift 
from the round-trip signal is given by 



%GM 



{{X^i^^+'sH^b^Sb-sH^l}']^..., (34) 



where the ellipses include terms of order {GMv /b){b/r) 
and higher order post-newtonian corrections. Note that 
the s''-' terms have vanished due to their oddness under 
Parity, just as they did for the time-delay formula. Also, 
the time dependence of is controlled by the behavior 
of the impact parameter vector b and its time derivative 

b. 



III. EXPERIMENTS 

In this section, the experimental implications of the 
results derived in Sees . Ill Cl and HTD] are examined in the 
context of key solar system experiments. We point out 



the basic features of the Lorentz-violating signals and 
contrast them with the GR case. Also, we describe how 
experiments could be used to probe various combinations 
of coefficients for Lorentz violation and estimate the level 
of sensitivity for each test. 



A. Preliminaries 

We work in a post-newtonian coordinate system that 
asymptotically coincides with the Sun-centered celestial- 
equatorial coordinate system adopted in most SME stud- 
ies [ll]. Space and time coordinates in this system are 
denoted with capital letters {T,X-') This approx- 

imation to an asymptotically inertial frame suffices for 
many SME experimental studies. Note that the Sun's 
center is in orbit around the barycenter of the solar sys- 
tem with a mean velocity about 1000 times smaller than 
the Earth's orbital velocity. Standard practice in solar 
system experiments is to adopt the Barycentric Celestial 
Reference System. For our purposes here in identifying 
the leading Lorentz-violating effects, it suffices to proceed 
in the Sun-centered frame and neglect the Sun's motion. 
However, in establishing beyond leading order corrections 
to the light travel time and Doppler observables in GR, 
the Sun's velocity can play a role [H, [l^l . 

To study the basic features of our results we focus 
on the solar conjunction time scale where the signals 
for Lorentz violation are near their maximum. In this 
scenario, where the light signal passes close to the Sun, 
we can assume approximately rectilinear motion for the 
Earth observer and the planet or spacecraft. The main 
changing variable in this case is the impact parameter 
vector [ssl [65| . Assuming rectilinear motion, we expand 
the impact parameter vector around its minimum value 
60 as 



5 = 60 + boT, 



(35) 



where 60 is the time derivative of the impact parameter 
vector evaluated at the conjunction time T = 0. Note 

that we also have b^ ■ b^ — Q. 

In many cases of interest, the time derivative of the 
impact parameter vector near the conjunction time can 
be written approximately as 



IpV -f leW 

R ' 



(36) 



where v is the Earth receiver's velocity and w is the ve- 
locity of the spacecraft or planet. All quantities on the 
right-hand side of Eq. ([36| can be determined from their 
definitions in Sec. Ill Bl and are evaluated at T = 0. Note 
that if the planet or spacecraft is many times further from 
the Sun than the Earth, so that R» le and L ~ i?, the 



primary contribution to ((36)) is from the Earth's velocity. 

The approximations described above will serve our 
purposes in exploring the features of the Lorentz- 
violating time-delay and Doppler signals. However, as 
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we discuss below in Sec. IIII C[ the more accurate results 
obtained in previous sections could be incorporated into 
a detailed computer code for a more rigorous approach. 
Furthermore, although wc focus below on the case where 
the central body is the Sun, many of our results can also 
be applied to the case where the Earth or other bodies 
produce the gravitational time delay and frequency shift. 



B. Time-delay and Doppler signals 

Adopting the solar system scenario described above 
where the central massive body is the Sun, we can es- 
tablish the general behavior of the time-delay formula. 
For definiteness, we adopt the post-newtonian coordi- 
nate gauge of Ref. setting ai = 02 = 1. Though 
this gauge differs from the standard harmonic gauge at 
0(3), for the 0(2) terms appearing in the time-delay ex- 
pression it is equivalent. Also, for times near conjunction, 
the gauge-dependent terms in ((26l) will be either approxi- 
mately constant or of order GMb/r or smaller, and hence 
negligible. 

The dominant contributions to the two-way time delay 
can be written as 



AT„ « 4GM 



(1 



In 



re+rp + R 
rr,+rr,-R 



(37) 

To illustrate the different functional dependencies of the 
terms in Eq. (j37[) we make use of the approximate ex- 
pression in ([55]) . Up to constants, the expression for the 
time delay becomes 



AGM 



TT 




+ S2- 



■bl + blT^ 



(38) 



where bo — \bo\ and 60 — \bo 
coefficients occurring in Eq. 



The two combinations of 
are given by 



Sl 



S2 



r'^ibfb^ 



-jku;k 



b^bo : 



(39) 



where 60 = bo /bo. 

There are three functions that appear in expression 
([55)) . The first term contains the standard logarithmic 
dependence present in GR, which is scaled by the rota- 
tional scalar combination of coefficients s'^'^ — + 



t;YY 



-zz 



s' ' + s'"'" . The second and third terms are controlled 
by the anisotropic combinations of coefficients si and S2 ■ 
To illustrate the typical behavior of the functions oc- 
curring in (|551). we plot them in Fig. [2] for the case of 
the Cassini experiment which took place near the solar 
conjunction on June 21, 2002. For this plot, we adopt 
the approximate values bo = l.Gi?©, 60 = 30km/s, and 
GM/c^ — 1.48 km, where Rq is the Sun's radius. 
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FIG. 2: The time-delay signals occurring in Eq. (|38|l near 
solar conjunction, plotted with the values for the Cassini ex- 
periment around June 21, 2002. The solid curve labeled GR 
gives the standard logarithmic dependence of GR, controlled 
by the combination 1 -I- s^"^. The curves labeled Ivl and lv2 
are the Lorentz-violating signals controlled by the combina- 
tions of coefficients si and S2, respectively. 



The logarithmic dependence of the time-delay signal 
is well known from GR [6^ . The two dashed curves in 
Fig. [2] represent departures from this standard behavior. 
In fact, part of the lv2 curve controlled by the combi- 
nation of coefhcients S2 produces an advancement of the 
light travel time rather than a delay. This may also occur 
with the lv2 curve controlled by the distinct combination 
of coefficients si , if the overall sign of this combination 
is negative. The peak values of the Ivl and lv2 curves 
are about 20/j,s in this example. Note that although we 
are effectively setting si = 1 and S2 = 1 for the pur- 
poses of plotting, no specific prediction is made here. As 
explained in the next subsection, these combinations of 
coefficients are expected to be much smaller than unity 
given current experimental constraints. 

It is also interesting to study the frequency shift that 
corresponds to the time-delay signal. For two combi- 
nations of coefficients for Lorentz violation in the grav- 
itational sector, this signal is enhanced over the time- 
delay signal. We examine the signal to 0(3) in the post- 
newtonian expansion and to leading order in b/r, assum- 
ing near-conjunction times. From Eq. (|34p the fractional 
frequency shift, expressed in the Sun-centered frame, is 
given by 



8GM 



[(1 



+ -s'H'b'')b 



-s'H'b''] 



(40) 



To see some of the features of the Lorentz-violating 
signals in (j40p . we use the approximate expression for 
the impact parameter vector psp . The expression for 
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the gravitational fractional frequency shift becomes 



— w 8GA/ 



-si 



(1 



blT 



blblT 



bl + blT^ 

bXihlT^ - hi) 



{hl + blT^Y 



S2- 



blT'^Y 



(41) 
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FIG. 3: The gravitational fractional frequency shift in Eq. 
(|40[) near solar conjunction, plotted with the values for the 
Cassini experiment around June 21, 2002. The solid curve 
labeled GR gives the standard 6/6 dependence of GR, con- 
trolled by the combination 1 + s^"^. The curves labeled Ivl 
and lv2 are the Lorentz-violating signals controlled by the 
combinations of coefficients si and S2, respectively. 



Just as in the time-delay case, three functions appear. 
We plot these functions in Fig. [31 again using values from 
the Cassini experiment and setting si = 1 and S2 = 1. 
The odd functional dependence of the signal controlled 
by the combination 1 + 's'^'^ is known [36, sl, The 
signal controlled by si resembles the GR case, though its 
peak size is reduced. The even functional dependence of 
the S2 signal is qualitatively different from the GR case. 
Note also that the maximum amplitude for this curve, 
which occurs at the conjunction time, is about twice that 
of the peak value for the GR curve. Also, as one can see 
qualitatively for each of the curves in Figs. [3| and [21 the 
Doppler signal is the neg ative of the time derivative of 
the time-delay signal [37[. 



Gravitational time-delay and Doppler tests provide no 
exception to this rule. In particular, the combinations of 
coefficients si and S2, controlling the Ivl and lv2 signals 
in Figs. [3 and [21 depend on the conjunction orientation 
of the experiment. To illustrate this, we include a sketch 
of the orientation of a typical experiment at the time 
of conjunction in Fig. [H This figure is oriented with 
the Sun-centered frame Z axis upwards, while re points 
in the ecliptic to the Earth's position. For experiments 
where the light signal comes within a few solar radii of the 
Sun, the spacecraft or planet position rp is only slightly 
inclined to the ecliptic. 




FIG. 4: Diagram illustrating the conjunction configuration 
of a typical time-delay or Doppler experiment in the solar 
system. The Sun-centered frame Z axis is shown along with 
the ecliptic plane (dashed line). 



As an example of this orientation dependence, con- 
sider the Cassini experiment in 2002. Near the time of 
conjunction, the Earth's velocity was pointing approx- 
imately along the Sun-centered frame X axis (i.e., the 
vernal equinox direction). Furthermore 6o was pointing 
very nearly perpendicular to the ecliptic. In this case the 
plane of the illustration in Fig. [H corresponds to the YZ 
plane with the X axis pointing out of the page. For this 
configuration we have 



(1,0,0), 
(0,0.4,-0.9). 



(42) 



C. Experimental analysis 

We discuss here key aspects of the experimental anal- 
ysis of the time-delay and Doppler signals for Lorentz 
violation in Eqs. (|37p and pp)) . Also, we make sensitiv- 
ity estimates for some key experiments. 

A ubiquitous feature of signals for Lorentz violation is 
the orientation dependence of observable signals [H, [l^l ■ 



This implies that the Cassini experiment is sensitive to 
the combinations of coefficients 

-s, « 0.25^^ + 0.85^^ - O.Ts^^ - 5^^, 

S2 ~ 0.4s^^ -0.9s-^^. (43) 

As another example, we suppose that the solar con- 
junction with the planet or spacecraft occurs near the 
vernal equinox. If this is the case, and the spacecraft or 
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planet is much further away from the Sun than the Earth 
and shghtly above the XY plane, we have 

^0 « (0,-1,0), 

bo « (0,0,1). (44) 

The Sun-centered frame coefficients for this scenario are 
given by 

52 « ~s^^. (45) 

Due to its scaling of the GR results in both the time- 
delay and Doppler signals, the rotational scalar combina- 
tion of coefficients s^^ is likely to be constrained at the 
same level as the PPN parameter 7, namely, parts in 10^. 
However, care is required since s'^'^ and 7 are not equiva- 
lent. In fact, the determination of the constant GM may 
correlate with s^^ . This is because s'^'^ occurs at 0(2) in 
Newtonian gravity [soj . For example, in orbital dynam- 
ics in the presence of s'^'' coefficients for Lorentz viola- 
tion, the basic Newtonian acceleration between two bod- 
ies is scaled by 1 -I- 5s^^ /3. If orbits are described as el- 
lipses with time-dependent orbital elements arising from 
perturbations to Newtonian gravity, the measured value 
{GM)jneas — n^a^, where n is the orbital frequency and 
a is the semimajor axis. Due to the presence of the s'^'^ 
coefficients {GM)^eas = GM{1 + 5s^^/3). We there- 
fore caution the reader that care is generally required in 
extracting constraints on s'^'^. 

To fit experimental data to the Lorentz- violating time- 
delay and Doppler signals, one could proceed by at least 
two methods. First, having already fit data to a GR 
signal, one could extract constraints on SME coefficients 
from time-delay or Doppler residuals. This could be ac- 
complished by using either the round-trip time-delay and 
Doppler formulas (|40|) and (jST]) or the less accurate ver- 
sions (|4ip and As mentioned before, one must bear 
in mind that the coefficients for Lorentz violation s'^'^ 
also affect orbital dynamics. The effects on the orbits 
of the planets due to the gravity-sector coefficients can 
be described as secular changes over orbital time scales, 
although oscillations can also occur [s^]. However, these 
effects could in principle be avoided with suitable filter- 
ing of the data if the focus is on the conjunction time 
scale b/v. 

Alternatively, detailed modeling of the time-delay sig- 
nal and the relevant orbital dynamics could be under- 
taken. In this case, the one-way formula in Eq. ([24|) . 
which is valid to 0{2) in the post-newtonian expansion 
and for times far from conjunction, could be used ap- 
propriately for both uplink and downlink. The full post- 
newtonian equations of motion for the Earth and space- 
craft or planet, and other relevant bodies that include 
the effects of the coefficients for Lorentz violation 's^'^ 
[sot , could be incorporated into the Orbital Determina- 
tion Program [6^. Indeed, data from past experiments 
using radar reflection from the inner planets [6^ [67| 



Experiment 


-TT 

S 


Si 


S2 


Ref. 


Time-delay signal 


Cassini 




10-^ 


10-* 


[38] 


Odyssey 




10-^* 


lo-''* 


[40] 


ASTROD 


io-» 


10-^* 


lo-'^* 


[42] 


BEACON 


10-9 


io-«* 


io-«* 




Doppler signal 


Cassini 


10-'^ 


10-* 


10-* 


[38] 


Odyssey 


10-^ 


10-^* 


IQ-'^* 


[40j 



TABLE I: Crude estimates of attainable sensitivities in some 
key experiments for the time-delay and Doppler signals. 



could be reanalyzed to search for SME coefficients via 
this second method described above. Although many of 
these past experiments lack data near conjunction, when 
the Lorentz-violating signals controlled by si and S2 are 
peaked, they could still be useful in measuring the ro- 
tational scalar combination s"^-^. Furthermore, detailed 
modeling may also reveal suppressed dependencies of the 
time-delay signal on combinations of s'^'' coefficients dis- 
tinct from s 1 and S2 ■ 

Regardless of the method adopted, we can make some 
reasonable estimates of the sensitivities achievable in ex- 
periments. We provide in Table [J estimated sensitivi- 
ties to the 3 dominant combinations of coefficients in 
the time-delay and Doppler experiments for some past 
and future experiments. We include the Cassini experi- 
ment and some key future tests. The estimates are order 
of magnitude only and are based on the peak values of 
the Lorentz-violating signals discussed above and the ap- 
proximate accuracy of each experiment referenced, when 
available. For example, the peak value of the S2 Doppler 
signal for the Cassini experiment is about IQ-^, while the 
Allan deviation for this experiment is about 10" |3§ |. 
indicating a sensitivity of parts in 10^. However, data 
from the time period when the S2 signal peaked in the 
2002 conjunction (T « in Fig. [3|) is not available, so 
the sensitivity to S2 is more likely to be parts in lO''. On 
the other hand, it appears likely that a suitable fitting 
of Cassini data could place the first constraints on the 
rotational scalar combination s-^-^ at the 10-^ level. For 
the time-delay signals, the sensitivity to the si and S2 
coefficients is reduced by about a factor of 10 or more, 
as indicated in Fig. [2l and this reduction in sensitivity is 
included in Table [J 

Proposals have been put forth for future experiments 
that measure to impressive accuracies the time delay 
from the Sun and even the Earth. We have included sen- 
sitivity estimates for the Odyssey, ASTROD, and BEA- 
CON experiments in Table |T1 Although in some cases 
it may be difficult to directly measure the fractional fre- 
quency shift (6^ . nonetheless we include some estimates 
in the table because of the possibility of increased sensi- 
tivity to SME coefficients from the Doppler signal over 
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the time-delay signal. The experiments in Table U are by 
no means an exhaustive list. Also of possible interest are 
proposals for measurin g th e light-bending effect such as 
SI M liil and LATOR {3^ , other proposed experiments 
[43I |47||. as well as existing accumulated data from Earth 
satellites [gO]- Though it lies beyond the scope of the 
present work, it would also be of interest to obtain the 
corresponding light-bending signal controlled by the s^" 
coefficients. 

Note that current constraints on the off diagonal com- 
ponents , , and s'^^ are at the level of 10~^ from 



atom interferometry S^] . Two combinations of these and 
other s"^^ coefficients are also constrained by lunar laser 
ranging at the 10^^*^ level [3l|. Thus, if future exper- 
iments can measure the peak behavior of the S2 set of 
coefficients in the Doppler signal to better than parts in 
10^, they may produce measurements of coefficients com- 
petitive with or better than previous experiments. The 
* label next to the estimated sensitivities in Tabic U in- 
dicates the requirement of measuring the peak behavior 
of the time-delay and Doppler signals. Finally we note 
that the S"^"^ coefficient does not appear at leading order 
in laboratory and orbital tests [30] and so time-delay and 
Doppler tests are likely to be among the most sensitive 
to this coefficient. 



IV. SUMMARY 

In this work, we have analyzed Lorentz-violating cor- 
rections to the gravitational time-delay and Doppler sig- 



nals in General Relativity, in the context of the gravita- 
tional sector of the minimal SME. We established general 
integral formulas for the deviation of a light ray from a 
straight line path that are valid in the linearized grav- 
ity limit. Our main results are analytical formulas for 
the light travel time and frequency shift for a light signal 
sent between two observers past a massive central body 
in the presence of gravity sector coefficients s^" . We ob- 
tained the one-way results in Eqs. ((24|) and ((32)l and the 
round-trip signals in Eqs. and (|34p . 



The Lorentz-violating signals were studied for solar 
system experiments involving light signals sent between 
the Earth and a planet or spacecraft near solar conjunc- 
tion. It was determined that the dominant signals are 
controlled by the combinations of coefRcients 1 -I- s"^"^, 
si, and S2- In terms of Sun-centered frame coefficients, 
the combinations si and S2 will vary for different experi- 
ments. We obtained sensitivity estimates for key existing 
and future experiments which are summarized in TablelH 
Time-delay and Doppler experiments could prove crucial 
in measuring the elusive scalar coefficient s'^'^ , to better 
than parts in 10^. Future highly sensitive time-delay and 
Doppler tests may be able to measure other coefficients in 
the subset s'^^ with sensitivities competitive with other 
existing experiments. 
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